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Abstract. Recently a method for adiabatic quantum computation has been 
proposed and there has been considerable speculation about its efficiency for 
NP-complete problems. Heuristic arguments in its favor are based on the un- 
proven assumption of an eigenvalue gap. We show that, even without the 
assumption of an eigenvalue gap, other standard arguments can be used to 
show that a large class of Hamiltonians proposed for adiabatic quantum com- 
putation have unique ground states. 

We also discuss some of the issues which arise in trying to analyze the 
behavior of the eigenvalue gap. In particular, we propose several mechanisms 
for modifying the final Hamiltonian to perform an adiabatic search with ef- 
ficiency comparable to that for 3-SAT. We also propose the use of randomly 
defined final Hamiltonians as a mechanism for analyzing the generic spectral 
behavior of the interpolating Hamiltonians associated with problems which 
lack sufficent structure to be amenable to efficient classical algorithms. 



1. Introduction 

Recently there has been considerable interest in a proposed scheme for adiabatic 
quantum computation ||^, |^ and speculation that it may even provide a mechanism 
for efficient solution of hard problems. Both the validity of the adiabatic theorem 
and the arguments |^ |^ for its efficiency depend on the existence of an eigenvalue 
gap. However, the existence of such a gap has not been proven. It has been 
conjectured |^ on the basis of numerical simulations and the so-called non-crossing 
rule. 

The adiabatic quantum algorithm is designed to take the ground state of an 
initial Hamiltonian Hq to that of a final Hamiltonian Hi using a linear interpolating 
Hamiltonian of the form H(s) = (1 — s)Hq + sHi. The quantum adiabatic theorem 



|14| , |19| , |20| , |24{ can be used to show that the efficiency of the adiabatic approxima- 
tion is 0(l/g^ij^) where gmin denotes the minimum energy gap g{s) = Ei{s) — Eq{s) 
between the ground and first excited states of H{s). Thus the efficiency of adia- 
batic quantum computation depends on how rapidly the eigenvalue gap decreases 
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as the size of the problem increases. The purpose of this note is not to resolve this 
question, but to discuss some of the issues that arise. 

Farhi, et al js], |^ (hereafter referred to as FGG) typically use an initial Hamil- 
tonian of the form 

(1.1) Ho = ^Y, ['^^(j) + = i XI ■ ■ ■ + ^] ■ ■ ■ ® ^ 

3 3 

with ttj a non-negative integer. The eigenstates of Hq are products of eigenstates 
of Gx- The so-called "computational basis" \ki . . . fc„) consists of products of eigen- 
states of the Pauli matrix a 2. (However, it is customary to identify each ki with an 
element of Z2 = {0, 1} rather than with the usual eigenvalues of ±1 or, equivalently, 
replacing cr^ by \[(Jz +-^] as has been done above.) In the computational basis, the 
eigenstates of Hq have the form 2^"/^ k \ki . . . kn) (with all signs +1 for 

the ground state). FGG define a non- negative final Hamiltonian Hi which is diag- 
onal in the computational basis and has the solution of some problem as its ground 
state, e.g., they encode the problem known as "3-SAT" in the computational basis. 
We note here only that this encoding results in a Hamiltonian of the form 

(1.2) ^1= E Ek,...kjki...k^)(ki...k^\ 

ki...kn 

where the ground state energy is zero and the other Ek-^...k„ are positive integers 
with an upper bound that is 0(n^). 

First, we point out that the non-crossing rule is completely inadequate for the 
purpose of establishing a gap. There are realistic models of physical systems which 
do exhibit crossings, despite the absence of a symmetry common to both the initial 
and final Hamiltonians Hq and Hi. The most well-known such system is H2 , the 
hydrogen molecule ion Jl^ . Another example is the Hubbard model for benzene 
|15{ . The non-crossing rule and the limitations on its applicability to adiabatic 
quantum computation are discussed in Appendix 

However, the most one could hope to gain from the non-crossing rule is unique- 
ness of the ground state. Fortunately, a standard argument based on the Perron- 
Frobenius theorem suffices for that purpose. This argument is sketched in 

the Section |^, but does not address the more fundamental issue of the size of the 
resulting eigenvalue gap. 

Some insight into the issues raised by adiabatic computation can be obtained 
by considering various strategies for replacing FGG's final Hamiltonian Hi by a 
modification whose ground state is the solution of Grover's search problem. We 
argue that if their algorithm is sufficiently robust to solve an NP problem in poly- 
nomial time 0(7VP), then a modification should be able to perform a successful 
search of an unordered list of M items in O {[log M]p) time, violating the conven- 
tional wisdom that Grover's algorithm (which requires 0(\/M) time) is optimal. 
A model problem suggested by FGG shows that this simplistic expectation need 
not hold; however, their example also suggests that gaps are associated with the 
presence of a symmetry, not its absence. 

The issues raised here involve questions in several subfields fields of physics 
and computer science. In order to make this note accessible to people with diverse 
backgrounds ranging from Schrodinger operator theory to computer science, two 
appendices are included — one on Grover's algorithm and one on the non-crossing 
rule. 
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2. Uniqueness of the Ground State 

Although proof of an eigenvalue gap is likely to be difficult, proving that systems 
of the type considered by FGG have a unique ground state is easier. It relies on 
a standard argument widely used to prove uniqueness (and positivity) of the 
ground state in a variety of systems, including quantum lattice models. There is 
no particular originality in the argument given below. We present it only in the 
hope of clarifying some issues, particularly the distinction between the uniqueness 
of the ground state and the existence of a lower bound on the size of the resulting 
eigenvalue gap. 



The ground state of the initial Hamiltonian (1.1) is easily seen to be unique and 
consists of products of ground states of ax{k). However, we wish to transform Hq 
to the computational basis of products \k1k2 ■ ■ -kn) of eigenstates of az{k). This 
is easily achieved using tensor products of the Hadamard transform. Instead of 
examining Hq itself, we consider the operator F — e^^° . The ground state of Hq 
is the eigenfunction corresponding to the largest eigenvalue of the matrix 



n 



(2.1) F = e-"° = g-afc[<T,(fc)+7]/2 ^Bi®B2...®Bn 

k=l 



where 



~ 2^ I 1 -1 I I e-«-/2 j I 1 -1 



(o o^ _ ^a,/2 f cosh^^ sinh^ 

^ ' \ sinh^f cosh^^ 

Now, since all elements of each Bk are strictly positive, their tensor product F = 
e~^° also has strictly positive elements. Hence, by the Perron-Frobenius theorem 
(Tt] , 21] the largest eigenvalue of F, and the ground state of Hq, is unique. 



We would like to know that the ground state remains unique for a Hamiltonian 
of the form H = {1 — s)Ho + sHi. For this, it suffices that Hi is diagonal (in the 
computational basis) and has a unique ground state. By the Lie- Trotter formula 



(2.3) e-" ^ lim (e'^"^e'^"'^ 

The effect of the diagonal matrix e~ is simply to multiply each row of e~~^^° 
by a positive number of the form e"'^^*/'". Hence the product e~"^^^e~~^^° also 
has positive elements and so does its m-th power. (Moreover, because e"'''^'*/™ 1 
as m — *■ 00, these positive elements do not become zero in the limit.) Thus, one can 
again apply the Perron-Frobenius theorem to conclude that the largest eigenvalue 
of e^^ is unique and, hence, H has a unique ground state if < s < 1- (The 
argument breaks down at s = 1 because e — I no longer has strictly positive 

elements off the diagonal and is completely decomposable.) 

Note that, even if Hi has a degenerate ground state, the interpolating Hamil- 
tonian iJ = (1 — s)Ho + tHi will still have a unique ground state for all < s < 1; 
however, the difference between the two lowest eigenvalues must ^ as s ^ 1. 
Thus, uniqueness of the ground state is a very different matter from an eigenvalue 
gap of minimal size. Indeed, the uniqueness argument above holds for models |Q, 
such as an adiabatic search for which the gap can be shown to decrease exponen- 
tially. 
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In addition to uniqueness, the Perron- Frobenius theorem imphes that the ground 
state has the form Ck\kik2 ■ ■ ■ kn) with strictly positive coefhcients > in the 
computational basis for all s € [0, 1). In the limit as s ^ 1 all but one of these 
Cfc 0. 

3. Adiabatic Search Algorithms 

As described in Appendix Grover's |T^ 13, 22 algorithm is designed to 
efficiently locate an unknown but identifiable target via the use of a unitary operator 
G which can be written as G = e"'^ where 



(3.1) A|fci...A:„)- 



\i{ki...kn) = {tl...tn) 

\ki...kn) otherwise 



One can easily implement an adiabatic search by choosing for the final Hamiltonian 
Hi in the FGG algorithm the Grover generator, A, above. The adiabatic evolution 
will take the ground state of Hq, namely, = 2""'/^ fc„ l^i ■ • • ^n) to the 

ground state of Hi — A, namely \ti . . However, FGG have shown that this 



process takes exponential time. The analysis can be simplified [25| by modifying 
the initial Hamiltonian Ho to reduce the analysis to a two-dimensional problem. 

The reduction to a two-dimensional problem, which plays a critical role in 
Grover's algorithm, is associated with a (2" — l)-fold degeneracy in the adiabatic 
Hamiltonian Hq and Hi. However, this is not at all essential for the success of an 
adiabatic search. All that is needed is that the ground state of Hi be the target 
state \ti . . . tn)- This suggests that one try to modify Hi so that its ground state is 
1^1 . . .t„), but the eigenvalue distribution of its excited states is similar to that of 
a final Hamiltonian known to have a gap. 

Suppose that a problem is known to have an efficient solution encoded in the 
final Hamiltonian Hi . Then setting H2 — GHi , yields a Hamiltonian identical to Hi 
except that the eigenvalue associated with the target state \ti. . . i„) is multiplied 
by —1. (Because G and Hi are both diagonal in the computational basis, they 
commute and H2 is self-adjoint.) Because Hi was defined to be non-negative, 
H2 has exactly one negative eigenvalue so that its ground state is now the target 
state \ti...tn)- Applying the adiabatic algorithm to the modified interpolating 
Hamiltonian i?2(s) = (1 — s)Hq + SH2 should take the ground state of Hq to the 
target state. Moreover, the only effect on the final Hamiltonian is to move one 
excited state below the previous ground state, without decreasing the final energy 
gap 5(1). A similar Hamiltonian which is non-negative could be constructed the 
replacement 

Eki...k„^ if fci...fc„ = ti...i„ 

Eki...k„ Eki...k„ + 1 otherwise. 



in (L2). In either case, one would not generally expect these modifications to 
significantly affect 5min(s), in which case the adiabatic search would be as efficient 
as the solution of the problem encoded in Hi . 

Thus, if Farhi, et al's projection of an efficiency of 0{nP) is correct for the 
problem encoded in Hi, then one would expect an efficiency of 0([log A^]^) for the 
adiabatic search of a list of = 2" items described above. However, this would 
imply that a quantum computer could search an unordered list in time 0([log A^]^), 
contradicting the conventional wisdom that a speed-up greater than 0{-\fN) is not 
possible 1^, ^, 22, 26]. This does not necessarily imply a contradiction. The 
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proofs that 0{\/ N) is optimal depend on assumptions about the nature of the 
"oracle query" used in the search. However, van Dam, Mosca and Vazirani | |25| 
have observed that the encoding of solutions of other problems, such as 3-SAT, in 
Hi implicitly assume the ability of the computer to perform more general queries. 
Thus, standard complexity query arguments can not rule out the possibility of 
polynomial time algorithms. 



After seeing a preliminary version of this manuscript, FGG |1C] pointed out 
that the change from Hi to H2 described above can have a dramatic effect on the 
gap. We describe their example in the next section. 

It should also be noted that our expectation that the change from Hi to H^ 
will not decrease the minimum gap is not based on the presumption that one is a 
small perturbation of the other. On the contrary, (as FGG | IC ] emphasized) the two 
Hamiltonians differ by a multiple of a projection, which can have a significant effect 
on the spectrum. Indeed, it is essential to our strategy to effect such a change in the 
final Hamiltonian. However, unless this also induces a change in the structure of 
the problem, such as a symmetry-breaking, this need not affect the generic behavior 
of the spectra of the interpolating Hamiltonian; in particular, it need not lead to 
an avoided crossing of the two lowest levels. 



4. A Separable Model 

Let Hi — \^j[(Tz{i) + I] and set Oj = 1 in H^. Then the interpolating 
Hamiltonian becomes 

(4.1) H{s) = Ho + s{Hi -Ho) = ^Yl [(^ " *)^-(^') + "^^^^'^ + ^ 

j 

This system is separable, and exactly solvable. Since [(1 — s)ax{j) + saz{j) + I 

has eigenvalues ^(1 ± y/1 — 2s + 2s^ ), the ground state of H{s) has energy 
^(1 — Vl ~ 2s + 2s2) and the gap g{s) — Vl — 2s + 2s'^ is independent of n with 
5min — The system also has a high level of symmetry, since H{s) commutes 
with elements of the symmetric group Sn- In fact, there is an additional accidental 

degeneracy so that the (fc + l)-st eigenvalue is k with a degeneracy of ^ ^ ^ for 

k = 0,l,2...n. 

Now if Hi is replaced by H2 = GHi, the symmetry is broken. However, there 
is still some symmetry and a high level of degeneracy. Essentially, only one state 
from each of the (n + 1) non-degenerate levels is affected and the problem reduces 
to an (71 + l)-dimensional one which can be analyzed explicitly and shown to have 
an exponentially decreasing gap. In effect, an eigenvalue can be associated with the 
target state and must "cross" the levels lying below it to reach the bottom of the 
spectrum when s = 1. (A more detailed examination shows that the levels become 
exponentially close and then bounce away with the target information transmitted 
to the lower level). 

Although this shows that the argument sketched in Section ^ above cannot be 
made rigorous, this model is not generic. The exponentially decreasing gap is the 
result of a symmetry breaking which should not occur when the Hamiltonian Hi 
has no symmetry to begin with. 
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5. Discussion 

Farhi, et al |^ analyze several other models for which the gap behavior can be 
calculated explicitly and shown to decrease slowly (i.e, polynomial in n). However, 
as they point out, these models all have a high level or symmetry or structure which 
would lead to efficient classical algorithms. In some cases symmetry allows a high 
level of degeneracy which permits one to squeeze 2" states into [0, 0{nP)] without 
forcing an exponentially decreasing gap. 

In adiabatic computation the typical choices for initial and final Hamiltonians 
have spectra with high degeneracies and consist of positive integers in a range that 
is polynomial in n. In general, the interpolating Hamiltonian H(s) breaks these 
degeneracies and must squeeze 2" distinct eigenvalues into a range of the form 
[0,O{nP)]. Thus, most of them must be exponentially close. 

It seem that a polynomial gap is more likely to be associated with the presence 
of symmetry, which allows high degeneracies, than with its absence. 

For excited states, it is irrelevant whether or not the observed mergings are 
"avoided crossings" or true crossings. The simultaneous coalescence of a large 
number of excited states is essential to the algorithm. The states of the initial 
and final Hamiltonians are both product states, those of the initial Hamiltonian 
are products of eigenstates of ax , while those of the final Hamiltonian are products 
of eigenstates of (and thus elements of the so-called "computational" basis). 
Therefore, the final ground state is always an evenly weighted superposition of all 
eigenstates of the initial Hamiltonian. Unlike standard applications of the adiabatic 
theorem, in which the main contribution comes from a few low-lying states, aZZof the 
excited states must contribute to the first-order correction. To do this, these higher 
excited states must get close in some sense. Fortunately, a quantum computer can 
make this first order correction efficiently, mixing in all 2" excited states, and this 
is where the method gets its potential power. However, in order that low order 
perturbation theory suffice, it is essential that the gap between the ground and first 
excited state not decrease too rapidly as n increases. 

The non-crossing rule, which is discussed in Appendix |^ is based on the belief 
that "accidents" are extremely rare so that such phenomena as persistent degener- 
acy, or level crossings do not occur without some underlying physical phenomenon 
(such as a symmetry) with implications for the associated mathematical model. 
This viewpoint would suggest that if the lowest gap is to decrease only polynomi- 
ally when the other eigenvalues are getting exponentially close, there must be some 
physical mechanism keeping them apart. We are skeptical that such a mechanism 
can be found for Hamiltonians which encode the solution of problems which do not 
have enough structure to yield efficient classical solutions. 

This raises another question. Is the spectrum of the interpolating Hamiltonian 
sensitive to the association of particular eigenvalues with particular eigenstates in 
the final Hamiltonian, or is it primarily dependent on the eigenvalue distribution? 
For Hamiltonians with a good deal of structure, the first situation clearly holds. 



and FGG have observed ||8|, 11| that the behavior of the eigenvalue gap can depend 
critically on the choice of initial Hamiltonian. However, one expects Hamiltonians 
which encode solutions to typical instances of classically intractable problems, to 
lack the structure needed for this sensitivity. 

Because of the difficulty in analyzing the behavior of the gap in problems with- 
out structure, it may be worth considering a randomly defined Hamiltonian, i.e.. 
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let the energy in ( pT^ ) have the form = /(fci . . . fc„) where / is a suitable 

random process. One can then ask if there is a sense in which the eigenvalue gap 
is "almost always" exponentially decreasing. If so, this would suggest that an ex- 
ponentially small gap is generic, and likely to occur for NP-complete problems. 
There is already an extensive literature [Q, ^ ^, on the spectral behavior of 
random Schrodinger operators, in which one obtains results about typical Hamilto- 
nians with certain properties rather than one with a fixed potential. Although the 
model Hamiltonians used here have quite a different structure, and may require the 
development of new techniques, this approach seems worth considering. 

Indeed, one could even define the final Hamiltonian for an adiabatic search by 
letting / be a random variable taking integer values in [1, n] and choosing 

Eki...k„ = ii ki . . . kn = ti . . .tn 

Eki...k„ = /(fci • ■ • kn) otherwise 



in (1.2) 
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Appendix A. Grover's Algorithm 



Grover's |22[ algorithm is designed to efficiently locate an unknown but 
identifiable target state \ti . . . tn). This state may be the key denoting the location 
of an item in an unsorted list (e.g., the analogue of the name in an alphabetized 
phone book associated with telephone number one has been given) or a state with 
a certain verifiable property, e.g., a representation of the factors of a given number 
or the solution of some NP-complete problem. 

Grover showed how to construct a unitary operator G whose effect is simply 
to multiply the unknown target state |ii . . . i„) by —1 and all others by +1. This 
operator can then be used to construct an algorithm which will find the target state 
with probability greater than i in 0{\/N), i.e., 0(2"/^) time for N — 2" states. 
The operator G is a unitary operator whose action on the computational basis is 
simply 

-\ki ...kn) if (fci . . . kn) = (ti... tn) 

\ki...kn) otherwise 



(A.l) G\k,...kn) = 



It is generally assumed that one has an "oracle" which can perform unitary oper- 
ations to determine whether or not a state has the desired property and outputs 
a function / whose value is 1 if the answer is yes and otherwise. This oracle is 
described by the unitary operator which takes 



(A.2) |fco) (»\ki... kn) ^ \ko ® f{ki . . . kn)) <»\ki... kn) 
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The action of this oracle when the first bit is in the Gx eigenstate 2 ^/^(|0) — |1) 
(for which we use the somewhat unconventional notation ll)^;) is then 

(A.3) |l),(^|fci...A:„) ^e"^«|l),(^|fci...fc„) 

which is exactly the effect of G when the ancillary initial bit ll)^; is omitted. The 
power of quantum computing is then exploited by applying G to superpositions of 
the form Cki...k„ \ki . . . kn) rather than to the individual states in the computa- 
tional basis. The analysis is facilitated by the realization that the problem can be 
reduced to a two-dimensional one in 



span||ii . . .i„), 



fei. 

for which |ti . . . t„) and |fci . . . fc„) form an orthonormal basis. 

^ ~ ^ fei...fc„#ti...t„ 

Appendix B. The Non-Crossing Non-Rule 

The so-called "non-crossing rule" is one of a number of physical principles which 
arise when symmetry ensures that a critical term in some expression, such as the 
leading term in a perturbation expansion, is zero. The remaining conditions needed 
to obtain a crossing, transition, etc. are then more easily satisfied. In the absence of 
the canonical conditions, such crossings and transitions are not truly "forbidden" , 
but are either rare events which result from an accidental confluence or the result of 
another physical circumstance which facilitates the satisfying of certain conditions. 

In the case of interest here, the canonical condition for crossing is that the 
Hamiltonians Hi and Hq both commute with the operators which generate a 
symmetry group Q (for which it is not necessary that Hi and Hq commute with 
each other). The irreducible representations of Q can then be used to classify the 
eigenspaces of Hi, Hq, and the interpolating Hamiltonian H{s) = (1 — s)Hi + sHi. 
The non-crossing rule asserts that one expects eigenvalues of II(s) to cross only if 
they belong to different irreducible representations. 

A similar situation arises in the Born-Oppenheimer approximation for the hy- 
drogen molecule ion iJ^ in which case the role of t is replaced by the internuclear 
distance R. Standard tables of atomic data, show many instances of crossings of 
states in the same symmetry class, in apparent violation of the non-crossing rule. 
In this case, the paradox was rather easily resolved by the discovery of a "hid- 
den symmetry" which can regarded as a generalization of that responsible for the 
well-known "accidental" degeneracies of states of different angular momentum (but 
same principle quantum number n) for hydrogen. However, the generator of this 



symmetry group is an operator (denoted F{R) by Judd |18[ ) which depends on the 
internuclear distance R. Although this gives a satisfactory physical explanation 
for the phenomena observed, it points out a difficulty with any attempt to make 
a mathematically precise theorem out of the "non-crossing rule" . If a crossing ex- 
ists, the degeneracy in H{s) always allows the formal construction of a suitable 



symmetry group |16| 



A less well-known example of violation of the non-crossing rule occurs in Heil- 



man and Lieb's study |15] of the Hubbard model for benzene. Moreover, in this 



case, a rather detailed analysis |lq] failed to locate a hidden symmetry, even one 



dependent on a parameter. In addition, their work found persistent degeneracies 
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unexplained by symmetry. This is more serious, as it is far less likely to be an 
artifact of the numerical methods used. It is worth quoting part of their discussion. 

The [non-crossing rule] depends [s] crucially on the interpretation of 
the word symmetry. The conventional meaning is that of a symmetry 
group independent of [a parameter] U; in this case the "proofs" are 
false. . . . [If] one allows symmetry groups that are [/-dependent, the 
"theorems" are mere tautologies, because . . . one can always invent, 
post hoc a [/-dependent group to account for any violations. 

One may ask what is wrong with the "proofs" quoted above. 
The fault lies not in the mathematics per se but in the assumptions 
used to connect the mathematics with the real world: First, in the 
natural sciences, two real numbers are never equal unless there is a 
physical reason for it; second, that reason must be the existence of a 
U-independent symmetry group. 

The first assumption has validity, but the second is merely a 
confession of ignorance ... 

Although the relevance of the non-crossing rule to adiabatic quantum computa- 
tion is questionable, the principles underlying it are not. If polynomially decreasing 
gaps are generic for systems in which exponentially many eigenvalues are squeezed 
into an interval that increases only polynomially, there must be a physical reason 
for this behavior. 
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